Aging and effective temperatures near a critical point 
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The orientation fluctuations of the director of a liquid crystal (LC) are measured after a quench 
near the Freedericksz transition, which is a second order transition driven by an electric field. We 
report experimental evidence that, because of the critical slowing down, the LC presents several 
properties of an aging system after the quench, such as power law scaling in times of correlation and 
response functions. During this slow relaxation, a well defined effective temperature, much larger 
than the heat bath temperature, can be measured using the fluctuation dissipation relation. 
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The characterization of the thermodynamic properties 
of out of equilibrium and of slow relaxing systems is 
an important problem of great current interest, which 
is studied both theoretically and experimentally. In this 
framework an important and general question, which has 
been analyzed only theoretically is the relaxation of a 
system, which is rapidly quenched exactly at the critical 
point of a second order phase transition [l], 0, S, 0, 0] • 
Because of the well known divergence of the relaxation 
time and of the correlation length the system presents a 
very rich dynamic, which has been named "aging at the 
critical point" . Indeed this relaxation dynamic presents 
several properties, which are reminiscent of those ob- 
served during the aging of spin glasses [Ell 11 II 0,0, [3 : 
the mean quantity has an algebraic decay and the corre- 
lation and response functions exhibit a power law scaling 
in time. Another important feature analyzed during the 
aging at a critical point concerns the equilibrium relation 
between response and correlation, i.e. the fluctuation 
dissipation theorem (FDT), which is not necessarily sat- 
isfied in an out of equilibrium system. This statement is 
relevant in the context of spin glass aging where a well de- 
fined (in the thermodynamic sense) effective temperature 
T e ff can be obtained using the so called fluctuation dis- 
sipation relation (FDR) . This relation leads to the intro- 
duction of the fluctuation dissipation ratio, X, between 
correlation and response functions @, 0] . Deviations of X 
from unity (X = 1 is the equilibrium value) may quan- 
tify the distance from equilibrium Equivalently an 
effective temperature T e s = T/X, larger than the heat 
bath temperature T, can be defined. Those ideas have 
been transferred to the aging at critical point, where for 
certain variables the FDR can be interpreted as a well 
defined effective temperature T e g [l|, 0] ■ 

The study of these analogies between the aging at the 
critical point and the aging of spin glasses is important 
because it allows to give new insight on the role of the 
quenched disorder on the above mentioned features and 
on the common mechanisms producing them. However 
these theoretical studies have been performed mainly on 
spin models and most importantly the reduced control 



parameter e has been set exactly equal to zero. Thus one 
is interested to know how general these predictions are 
and also one may wonder whether those predictions can 
be observed in an experimental system where the exact 
condition e = can never be reached. 

The purpose of this letter is to analyze, within this 
theoretical framework, an experiment on the relaxation 
dynamic close to the critical point of a liquid crystal in- 
stability, which is, at a first approximation, described by 
a Ginzburg-Landau equation and has the main features 
of a second order phase transition. The main result of 
our investigation is that after a quench close to the crit- 
ical point the effective temperature T e g, measured from 
FDR, has a well defined value larger than the thermal 
bath temperature. 

This kind of experimental test is important because 
the properties of X and T c g have attracted much interest, 
since they suggest that a generalized statistical mechanics 
can be defined for a broad class of non equilibrium phe- 
nomena. Several theoretical models or numerical simula- 
tions show such a behavior . However in experiments 
the results are less clear and the relevance for real materi- 
als of the T e g as defined by FDR is still an open question, 
which merits investigation [ljj, 11, E3, EH- For exam- 
ple, there are important differences, which are not un- 
derstood, between supercooled fluids polymers 12 1, 



gels [111] and spin glasses [13J. The point is that on this 
subject it is very difficult to find simple theoretical mod- 
els, which can be directly compared with experiments. 
Thus the experimental study of the relaxation dynamics 
close to the critical point is very useful in this sense. 

The system of our interest is the Freedericksz tran- 
sition of a liquid crystal (LC), subjected to an electric 
field E 14, HI]. In this system, we measure the variable 
C, which is the spatially averaged alignment of the LC 
molecules, whose local direction of alignment is defined 
by the unit pseudo vector n. Let us first recall the gen- 
eral properties of the Freedericksz transition. The system 
under consideration is a LC confined between two par- 
allel glass plates at a distance L = 9 ^m. The inner 
surfaces of the confining plates have transparent Indium- 
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Tin-Oxyde (ITO) electrodes, used to apply the electric 
field. Furthermore the plate surfaces, are coated by a 
thin layer of polymer mechanically rubbed in one direc- 
tion. This surface treatment causes the alignment of the 
LC molecules in a unique direction parallel to the surface 
(planar alignment),! e. all the molecules have the same 
director parallel to the z-axisflij]. The cell is next filled 
by a LC having a positive dielectric anisotropy (p-pentyl- 
cyanobiphenyl, 5CB, produced by Merck). The LC is 
subjected to an electric field perpendicular to the plates, 
by applying a voltage V between the ITO electrodes, i.e. 
E = V/L. To avoid the electrical polarization of the 
LC, wc apply an AC vol tag e at a frequency fy = 1 kHz 
(V = V2V cos(2ir f v t)) [3 HE]. More details on the ex- 
perimental set-up can be found in ref . [3, EH ■ When Vq 
exceeds a threshold value V c the planar state becomes 
unstable and the LC molecules, except those anchored to 
the glass surfaces, try to align parallel to the field, i.e. the 
director, away from the confining plates, acquires a com- 
ponent parallel to the applied electric field (z-axis). This 
is the Freedericksz transition whose properties are those 
of a second order phase transition [lj, ll_5[ . For Vq close 
to V c the motion of the director is characterized by its an- 
gular displacement 9 in a;z-plane (9 is the angle between 
the x axis and ft), whose space-time de pen dence has the 
following form : 9 = 6 (x, y, t) sin (f) [3, 0, E$k- If #0 
remains small then its dynamics, neglecting the (x, y) de- 
pendence of 9q, is described by a Ginzburg-Landau equa- 
tion : 
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(n + e + i)el + rj 



(1) 



where e = p%- 



1 is the reduced control parameter. The 
characteristic time To and constant k depend on the LC 
material properties, rj is a thermal noise delta-correlated 
in time [l9(. 

We define the variable C as the spatially averaged align- 
ment of the LC molecules and more precisely : 

q = Ta\\ dxdy l (1_ri ' )42 (2) 

where A = ttDq/A is the area, in the (a;, y) plane, of the 
measuring region of diameter D , which is about 2 mm 
in our case. If 9 remains small, £ takes a simple form in 
terms of 9q, i.e. n x = cos(9) and £ = 1/A JJ A dxdy9'Q. 
It has been shown that £ is characterized by a mean value 
(C) oc e, by a divergent relaxation time r = T /e an d by 
fluctuations, which have a Lorcntzian spectrum 18l. |21|. 
The measure of the variable £ relies upon the anisotropic 
properties of the LC, i.e. the cell is a birefringent plate 
whose local optical axis is parallel to n. This optical 
anisotropy can be prec isely measured using a polariza- 
tion interferometer [201 ] , which has a signal to noise ratio 
larger than 1000 (see ref. [3, 13 for details). 

In this letter, we consider the dynamics of ( as a 
function of time after a quench from e = e\ > 0.3 to 
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FIG. 1: a) Mean relaxation of the system, log(£(t)), as a 
function of time after a quench (o) from e\ = 0.3 to eo = 0.01. 
The continuous line represents the theoretical response based 
on eq. [T]with (£) — ipo- b) Mean relaxation of the system as 
a function of the reduced time (t + r m ). in both figures the 
dashed lines represents the theoretical relaxation at e = 0. 



e = eo << 0.1. The fact that the control parameter is the 
electric field allows one to have extremely fast quenches, 
typically 1 ms. As the typical relaxation time r of ( 
is about 10s, at eo — 0, this means that we can fol- 
low the out of equilibrium dynamics for about four or- 
ders of magnitude in time, which is really comparable to 
what is done in real glass experiments after a tempera- 
ture quench, where the dynamics is usually followed for 
about 3 or 4 decades in time. The advantage here is that 
this relaxation dynamics, which lasts only a few seconds, 
allows us to repeat the experiment several times and to 
perform an ensemble average (indicated by < • >) of the 
measured quantities. We consider first the specific case 
of a quench from ei = 0.3 to eo = 0.01. The typical 
mean value < £(i) > of £ after this quench is plotted in 
fig. [I] as a function of time t (t = is the time when the 
quench has been performed). This mean dynamics of £ is 
obtained by repeating the same quench 7000 times. The 
behavior of < £ > remains constant for a certain time and 
then slowly relaxes (see fig. QJi). Above a characteristic 
time, which is about 0.22 s in our case, the relaxation 
becomes exponential. 

To understand this behavior, we decompose the dy- 
namics of #o in a mean dynamics after the quench and its 
fluctuations, i.e. 9q = ipo(t) +59o, where xpo(t) = (9o(t)). 
The dynamics of "fo obtained from the analytical solu- 
tion of cqfTJis shown in fig[T]as a continuous line, which 
perfectly agrees with the measured dynamics. The ex- 
perimental data and the solution of cqQ] have two well 
distinguish limits : for t 3> r = To/(2e), the relax- 
ation is exponential (r ~ 0.22s at eo = 0.01) ; for 
Kte To/(2eo), the dynamics of ipo is almost algebraic 

-1/2 



(t + T m ) 



with 



(k+1)>0o(O) 2 - 

This behavior, plotted in fig. [T] as a dashed line, is iden- 
tical to the relaxation at the critical point [J]. Thus 
the system should present aging phenomena in the range 
t < r, which is the interval where the dashed line follows 
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FIG. 2: a) Correlations functions C$(t,t w ) as a function of 
t - t w for different fixed t = 0.06 s (o), 0.2 s (□), 0.26 s (o), 
0.28 s (x) and 0.31 s (+) after the same quench of fig.Q] Inset: 
The correlation functions have a simple master curve obtained 
by plotting Cg(t,t w )/Ce(t,t w ) versus (t—t w )/(t + T m ). Notice 
that (t + T m ) is the reduced time used in fig[T] 
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FIG. 3: a) Parametric plot of integrated response versus cor- 
relation for t = 0.4 s (o), 0.5 s (□), 0.6 s (o) and 0.7 s (+). 
Continuous line are linear fits. Dark line represents the FDT 
at equilibrium. Inset : the characteristic time t* is plotted as 
a function of t. 



the experimental data in fig. [T] 

In order to measure T c g as denned by the FDR, we 
need to measure both the correlation function Cg(t, t w ) = 
( S6o(t) 58o(t w )} of the thermal fluctuations 89q of 9q 
and the response function x{t,t w ) °f at time t to a 
perturbation, given at time t w , to its conjugated variable 
F, with t w < t. More precisely, from ref . [6|, IS)] the effective 
temperature is given by: 
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(3) 



where ks is the Boltzmann constant, T the heat bath tem- 
perature, i.e. T G ff = T/X(t,t w ). 

To compute x(t->t w ) and Cg(t,t w ) in our experiment 
one has to consider that the measured variable is £. 
The relationship between the fluctuations of £ and 9q 
has been discussed in details in ref.[l7j, thus we recall 
here only the useful results. As the area of the mea- 
suring beam is much larger than the correlation length, 
the global variable measured by the interferometer is 
C = 4 // O^dxdy ~ ifjQ+2tlj S9 . Thus the mean value of 
£ is (C(t)) = ipo(t) and the fluctuations of S( of C can be 
related to the fluctuations of 8q : 5((t) = 2ipo(t)S9o(t). 

The autocorrelation function of 9o is obtained us- 
ing the values of ipo(t) and tpo(t w ), i.e. Cg(t,t w ) = 
{((t)((t w ))/(4Mt)Mtw))- The Cg(t,t w ), measured at 
various fixed t, are plotted as function of t — t w in figiU 
We see that a simple scaling of the correlation func- 
tion can be obtained by plotting them as a function of 
(t — t w )/(t + r m ). Thus in agreement with theoretical 
predictions, our system presents, during its slow relax- 
ation after the quench close to the critical point, algebraic 
decay of mean quantities and the scaling of correlation 
functions, which are very similar to the main features of 
spin glass aging (see for example @, 0] ) . 

The response function is obtained by applying a small 
change of the voltage Vb, which modifies the control pa- 
rameter, i.e. e = eo + <5e. In ref.fl?]] we have shown that 
the external torque T ext , i.e. the conjugated variable of 



#o, associated to 5e is equal to 

r GXt = 2B Se Mt) (l 
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We separate 8o(t) into the average part ifjo(t), solution 
of eq. [U with 77 = 0, and a deviation A(t) due to 
Se : 9o(t) = ipo{t) + A(t). We define the linear re- 
sponse function R{t, t w ) =< A > /r cxt of 9q to T ext us- 
ing a Dirac delta function for 8e at an instant t w , i.e. 
Se(t) oc 5(t — t w ). In the experiment the Dirac delta 
function is approximated by a triangular function of am- 
plitude 5eo and duration r r ~ 2ms << t, specifically: 
Se = 5e a (l - 2|t - t w \/r r for \t - t w \ < r r /2 and Se = 
for \t — t w \ > r r /2, with T r /r ~ 10 -4 . The measured 
quantity is the response R<^,Se =< A( > /Se, of ((t) to 
Se(t w ). As AC(t) = 2if) (i)A it follows that the linear 
response of 9q to T ext is: 



R(t, t w ) 
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where B = Air 2 Ki/(4:L), and K\ a LC elastic con- 
stant. Thus inserting in eqf5] the measured values of 
R^ t Se(t,tw), of ipo(t) an d of i/jo(t w ) we can measure 
R(t,t w ) and by numerical integration we finally obtain 
X(t,t w ) = f t w R(t,t')dt' '. The accuracy of this procedure 
has been checked at equilibrium in ref. [171 ] , where we have 
shown that the measured Cg a and x(*>^u) verify FDT. 

The integrated response is related to the correlation by 
the FDR relations (see eq.[3]). Due to the non-equilibrium 
process, it is not equivalent whether the parameter is the 
waiting time t w or the observation time t (see ref.Q for 
a detailed discussion on this point). Applying the cor- 
rect procedure 0, @| , we keep the time t constant and we 
vary t w , i.e. < t w < t, and we repeat the procedure at 
various t. Following ref.0], we study the FDR during the 
relaxation process by plotting the integrated response x 
as a function of Cg. These FDR plots can be seen in 
fig. [3] for four characteristic values of t. When the sys- 
tem is in equilibrium the FDR plot is a straight line with 
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slope — l/(fcsT) (continuous line in fig. [3]). The curves in 
fig. [3] for each t are composed by two straight lines and 
are remarkably similar to those predicted in theoretical 
models Q and seen in experiments [f3| of spin glasses 
aging. When t w is close to t, i.e. large values of C(t, t w ), 
the FDT is satisfied, i.e. the experimental points are on 
the continuous line in fig. [3] For C(t,t w ) smaller than a 
value C* , which depends on t, the FDT is not satisfied. 
However x{t, t w ) remains linear in C(t, t w ) but the slope, 
which does not depend on t, is much smaller than its 
value at equilibrium. This means that the slow modes of 
the systems have a very well defined T e ff, independent 
of t, higher than the temperature of the bath. Precisely 
we find, for C(t,t w ) < C* , X = T/T cS = 0.31, which is 
remarkably close to the asymptotic value X = 0.305 an- 
alytically Q and numerically Q estimated for the small 
wave vector modes of an Ising model quenched at the 
critical point using a different quenching procedure. This 
analogy is important because our measurement, being an 
integral over the measuring volume, is more sensitive to 
the fluctuations of the long wavelength modes. In fig[3] 
we see that C* is a decreasing function of t. Let us de- 
fine t* as the characteristic time associated with C* . In 
the inset of fig. [3j we have plotted the value of t* as a 
function of t. We find that t* ~ for t < r = 0.22s. 
This indicates that, in agreement with refT0, 0, HI, the 
violation occurs for all times at 6q = 0. In our case it is 
cut by the exponential relaxation, which starts at t ~ r. 
For t > r, t* = t — t is linear in t, which also agrees 
with the theoretical picture of the FDR Indeed t* 
separates the equilibrium part and the aging one and the 
ratio t*/t = 1 — r/i, for t > r, defines how large is the 
equilibrium interval with respect to the total time. No- 
tice that at eo = 0, r = oo and the equilibrium interval 
does not exist Thus the values of eo and e\ affect 
only the amplitude of the region where the out of equi- 
librium is observed but not the value of X 22j . We also 
want to stress that this very clear scenario, with a well 
defined T e g independent of t, is obtained only if the right 
procedure with fixed t is used [D, 0] . The results are com- 
pletely different if t w is kept fixed. In such a case T e ff is 
a decreasing function of t w [22j • The non-commutability 
of the two procedures, i.e. either with t or with t w fixed, 
shows the non ergodicity of the phenomenon and opens a 
wide discussion on how to analyze the data in more com- 
plex situations, such as those where T e ff is measured in 
real materials and the procedure with t w fixed is often 
used [HI [II [II, [HI . 

In conclusion we have shown that a LC quenched 
close to the critical point of the Freedericksz transition 
presents aging features, such as a power law scaling of 
correlation functions and the appearance of a well de- 
fined T e ff > T. Most importantly we show the existence 
of a well defined T c g if the the right procedure with fixed 
t is used. What is very interesting here is that although 
we are not exactly at eo = we observe a large interval 



of time,0 < t < t where the predicted aging at critical 
point can be observed. The results plotted in figl3] agree 
with those of ref. 0, Hi but are different from those pre- 
dicted for mean-field 4|. This opens the discussion for 
further theoretical and experimental developments. It 
also shows that the study of the quench at critical point 
is an interesting and not completely understood problem 
by itself 0,0]. 
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